Testing ACDM with the Growth Function 5(a): Current Constraints 
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We have compiled a dataset consisting of 22 datapoints at a redshift range (0.15,3.8) which can be 
used to constrain the linear perturbation growth rate / = jfr^ . Five of these data-points constrain 
directly the growth rate / through either redshift distortions or change of the power spectrum with 
redshift. The rest of the datapoints constrain / indirectly through the rms mass fluctuation as{z) 
inferred from Ly-a at various redshifts. Our analysis tests the consistency of the ACDM model 
and leads to a constraint of the Wang-Steinhardt growth index 7 (defined from / = f2 m (a) 7 ) as 
7 = 0.67lo J?- This result is clearly consistent at la with the value 7 = yy = 0.545 predicted 
by ACDM. We also apply our analysis on a new null test of ACDM which is similar to the one 
recently proposed by Chiba and Nakamura (arXiv:0708.3877) but does not involve derivatives of 
the expansion rate H(z). This also leads to the fact that ACDM provides an excellent fit to the 
current linear growth data. 

PACS numbers: 



I. INTRODUCTION 

Most of the current observational evidence for the ac- 
celerating expansion of the universe and the existence of 
dark energy comes from geometrical tests that measure 
directly the integral of the expansion rate of the universe 
H(z) at various redshifts. Such tests include measure- 
ments of the luminosity distance by using standard can- 
dles like type la supernovae (Snla) [l| and measurement 
of the angular luminosity distance using standard rulers 
(last scattering horizon scale [2| , baryon acoustic oscilla- 
tions peak Q etc). Even though these tests are presently 
the most accurate probe of dark energy, the mere deter- 
mination of the expansion rate H(z) is not able to pro- 
vide significant insight into the properties of dark energy 
and distinguish it from models that attribute the accel- 
erating expansion to modifications of general relativity. 
The additional observational input that is required is the 
growth function 6(z) = -^(z) of the linear matter den- 
sity contrast as a function of redshift. The combination 
of the observed functions H(z) and 5{z ) can provide sig- 
nificant insight into the properties of dark energy (e.g. 
sound speed, existence of anisotropic stress etc) or even 
distinguish it from modified gravity theories 

The observational description of the expansion rate 
H(z) is usually made through the use of parameteriza- 
tions of the effective equation of state of dark energy 
Pdcl ) z \ which is related to H(z) as 
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The most commonly used such parametrization is the 
CPL parametrization 0, 0] 



w(z) = wq + w\ 



l + z 



(1.2) 



The cosmological constant (ACDM ) corresponds to pa- 
rameter values wo — — 1, w\ — and is not only con- 
sistent with all current geometrical tests but it is also 



favored 0, 0] by most of them compared to other param- 
eter values. 

The corresponding parametrization of the linear 
growth function S(z) can be made efficiently by intro- 
ducing a growth index 7 defined by 



din 8(a) 
din a 



where a 



— is the scale factor and 
1+2 



t (a) = 



H 2 (a) 
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This parametrization was originally introduced by Wang 
and Steinhardt [8| (see also [9( for more recent discus- 
sions) and was shown to provide an excellent fit to 
corresponding to various general relativistic cosmological 
models for specific values of 7. In particular, it was shown 
that for dark energy with slowly varying w(z) ~ wq the 
parameter 7 in a flat universe is 



7: 
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which for the ACDM case (w = —1) reduces to 7 = yy. 
It is therefore clear that the observational determination 
of the growth index 7 can be used to test ACDM . 

The observational determination of the growth index 7 
requires knowledge not only of S(z) but also of H(z) and 
ftom ( see equation (|1.4| ). It is possible however to con- 
struct more direct tests of ACDM which require knowl- 
edge only of 6(z) and H(z). Such a null test of ACDM 
was recently proposed in Ref. [1 01 ] where it was suggested 
that the validity of ACDM requires the following equality 
of observables 



(H(z)yHjy 

(l + z) 2 ^'(0) 2 
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where ' = 4^. In fact, as discussed section III, there is an 
improved version of this null test that does not involve 



2 



derivative of the observable H(z) and therefore it is less 
prone to observational errors. 

Both of the ACDM tests discussed above (the growth 
index 7 and the null test) require observational deter- 
mination of 8{z). There are several observational ap- 
proaches that can lead to the determination of 5(z). For 
example, redshift distortions of galaxy power spectra [ill ] , 
the rms mass fluctuation as(z) inferred from galaxy and 
Ly — a surveys at various redshifts [Hj]-[l3|], weak lens- 
ing statistics 14 1, baryon acoustic oscillations [ll|, X- 
ray luminous galaxy clusters [l|| , Integrated Sachs- Wolfe 
(ISW) effect [13] etc. Unfortunately, the currently avail- 
able data are limited in number and accuracy and come 
mainly from the first two categories. They involve signif- 
icant error bars and non-trivial assumptions that hinder 
a reliable determination of 5{z). In addition, a large part 
of the available data are at high redshifts (z > 1) where 
ACDM is degenerate with most other dark energy mod- 
els since dark energy is subdominant compared to matter 
at high redshifts in most models. 

Nevertheless, it is still instructive to consider the 
presently available data to investigate the possible weak 
constraints that can be imposed on ACDM . Such a task 
serves two purposes 

1. It can be used as a paradigm for the time when 
more accurate data will be available 

2. It can provide constraints from a dynamical test 
which are orthogonal and completely independent 
from the usual geometrical tests. 

Thus, in what follows we use a wide range of presently 
available data on both redshift distortions and o%(z) to 
determine the observed growth index 7 and test ACDM 
in two ways 

1. Comparing the measured value of 7 with the 
ACDM prediction 7 = £ 

2. Implementing a new null test that exploits the con- 
sistency between H(z) and S(z) in the context of 
ACDM . 

The structure of this paper is the following: In the next 
section we present the dataset we have compiled and use 
it to determine the best fit value of 7 under the assump- 
tion of a ACDM background. In section III we derive 
a new null test for ACDM and apply it using the fit 
performed in section II. Finally, in section IV we con- 
clude, summarize and present the future prospects of the 
present study. 



II. FITTING THE GROWTH INDEX 

According to general relativity, the equations that de- 
termine the evolution of the density contrast 6 in a flat 
background consisting of matter with density p m and 



dark energy with pd e 



I'dc 



are of the form 



S + 2-6 = 4irGp m 5 
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It is straightforward to change variables in eq. (|2.ip from 
t to In a = H-^) to obtain § 



(ln5)"+(ln<5)' 2 +(m<5)' 
where we used (I2.3|) and 
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as in (|1.4p . A further change of variables from In a to 
^m(a) can be made by considering the differential of eq. 
(|2.5p and using energy conservation (dp = — 3 (p+p)d In a) 
leads to 



df£ m = 3wfl m (l — Q m )dliia 
Using (HU) in ([21| we get @ 
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where we have set 
Using the ansatz 
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in eq. (|2.7p an d expanding around Q m — 1 (good ap- 
proximation especially at z ~ 1) we find to lowest order 

_ 3(w - 1) 



6w — 5 



(2.10) 



which reduces to 7 = yj for ACDM (w = —1). 

Equations (|2.9p and (|2.10p provide excellent approxi- 
mations to the numerically obtained form of f(z). This 
is demonstrated in Fig. 1 where we plot the numerically 
obtained solution of eq. (|2.4|1 for the normalized growth 



( \ - S ^ 



(2.11) 



in the case of ACDM (Qom = 0.3) along with the corre- 
sponding approximate result 



(2.12) 
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Numerical Solution 



Approx. Solution (y= — ) 



TABLE I: The currently available data for the parameters /3 
and b at various redshifts along with the inferred growth rates 
and references. Notice that Ref. [23|] only reports the growth 
rate and not the /3 and b parameters since the growth rate was 
obtained directly from the change of power spectrum Ly — a 
forest data in SDSS at various redshift slices. 



z 


P 


b 


fobs Ref. 


0.15 
0.35 
0.55 
1.4 
3.0 


0.49 ±0.09 
0.31 ±0.04 
0.45 ± 0.05 

o.6o±°:li 


1.04 ±0.11 
2.25 ±0.08 
1.66 ±0.35 

1.5 ±0.20 


0.51 ±0.11 [11], [24] 
0.70 ±0.18 [25] 
0.75 ±0.18 [26] 
0.90 ±0.24 [27] 
1.46 ± 0.29 [23] 



FIG. 1: The numerically obtained solution of eq. (|2.4|l for 
the normalized growth of eq. (|2.11|l in the case of ACDM 
(Oom = 0.3) (black dashed line) along with the correspond- 
ing approximate result with j = A obtained from eq. (|2.13[) 
(blue continuous line). The agreement between the two ap- 
proaches is excellent. 



(for a pedagogical discussion see [20]). The parameter 
is related to the growth rate / as 
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with 7 = ji obtained from 

fin \ dS / da n ( ^7 



(2.13) 



The difference between the two approaches is less than 
0.1%. 

Despite of the impressive agreement between numeri- 
cal result and analytical approximation, there has been 
a recent attempt [18j to improve further on the analyt- 
ical approximation by considering an expansion of 7 in 
redshift space up to first order in z 



7 = 7o ± 7o z 



(2.14) 



where 70 and 7 are constants. Using (|2.14j) in (|2.9j) and 
(|2.7p it may be shown that [18[ 
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2 m + 2 



0.3 becomes 



which for 70 = yr, iwrj ■ 
7q = —0.012. Since j' is very small for ACDM we set it 
to zero in what follows and assume 7 = 70 = constant. In 
the Appendix however, we generalize our fit using equa- 
tion (|2.14[) and show that the introduction of the new 
parameter "f' Q increases significantly the errors and the 
allowed parameter region at la. 

Our goal in this section is to fit the parameter 7 using 
observational data and compare it with the value 7 = 
■A- of ACDM . The most useful currently available data 
that can be used to constrain 5(z ) (and 7) involve the 
redshift distortion parameter j3 [191 ] observed through the 
anisotropic pattern of galactic redshifts on cluster scales 



where b is the bias factor connecting total matter pertur- 
bations S and galaxy perturbations 5 g (b = ^-). 

The currently available data for the parameters /3 and 
b at various redshifts are shown in Table I along with the 
inferred growth rates and references. This is an expanded 
version of the dataset used in Ref. [2l[ where a similar 
analysis was performed using a different parametrization 
suitable for modified gravity models. Notice that the 
Refs. of Table I have assumed ACDM (with fl 0m = 0.3) 
when converting redshifts to distances for the power spec- 
tra and therefore their use to test models different from 
ACDM may not be reliable. In addition, as pointed out 
in Re f. [22l| , the points of Table I obtained from Refs [25| , 
(26|, |27j correspond to measurements of (3(z) but the es- 
timate of the bias b is made by comparing numerically 
simulated power spectra with the observed ones. Since 
the numerically simulated power spectra have assumed 
a ACDM cosmology, the resulting f b s = (i b should be 
interpreted carefully can only be used to test the consis- 
tency of ACDM cosmology. Given this ambiguity we have 
evaluated the best fit value of the index 7 both with and 
without the three points discussed above (see equations 
(|2~23D and 

In the same Table we also show the growth rate ob- 
tained in Ref. [23[ which does not rely on j3 but is based 
on a different strategy, namely by finding directly the 
change of power spectrum Ly — a forest data in SDSS 
at various redshift slices. This point (in addition to the 
first one) has been used previously by other authors (see 
Ref. [21]) in a similar context. 

Using the data of Table I we can perform a maximum 
likelihood analysis in order to find 7 and check its con- 
sistency with the ACDM value yj. We thus construct 



X/(^0m,7) = X! 



fobs(zi) - fth{Zi,l) 



°7o 



(2.17) 
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where f b s and a fobs are obtained from Table I while 
fth{ z iil) is obtained from eq. (|2.9[) . 

An alternative observational probe of the growth func- 
tion 5{z) is the redshift dependence of the rms mass fluc- 
tuation erg (z) defined by 



a 2 (R,z) = / W 2 (kR)A 2 (k,z)^ 



with 



W(kR) 
A 2 (kz) 



sin(fc-R) cos(fci?)\ 
(fci?) 3 (kR) 2 ) 
4:TTk 3 Ps(k,z) 



(2.18) 

(2.19) 
(2.20) 



with R = 8h~ 1 Mpc and Ps(k,z) the mass power spec- 
trum at redshift z. The function cr s (z) is connected with 

S(z) as 



which implies 



0) 



(2.21) 



sth (21,22) 



a 8 (z 2 ) S(z 2 ) 



(2.22) 



TABLE II: The currently available data for the rms fluctua- 
tion ag(z) at various redshifts and references. Notice that the 
data from Ref. [2^] were obtained using the normalized by a 
growth factor 8 and in our analysis we took this into account. 



z 


<X8 


0" CT8 


Ref. 


2.125 


0.95 


0.17 


[12] 


2.72 


0.92 


0.17 




2.2 


0.92 


0.16 


[13J 


2.4 


0.89 


0.11 




2.6 


0.98 


0.13 




2.8 


1.02 


0.09 




3.0 


0.94 


0.08 




3.2 


0.88 


0.09 




3.4 


0.87 


0.12 




3.6 


0.95 


0.16 




3.8 


0.90 


0.17 




0.35 


0.55 


0.10 


[28] 


0.6 


0.62 


0.12 




0.8 


0.71 


0.11 




1.0 


0.69 


0.14 




1.2 


0.75 


0.14 




1.65 


0.92 


0.20 





where we made use of eq. (|2.12|) . Most of the currently 
available datapoints ag(zi) originate from the observed 
redshift evolution of the flux power spectrum of Ly — a 
forest [HI, EH, HH . These datapoints are shown in Table If 
along with the corresponding reference sources. Notice 
that the data from Ref. [28| were obtained using the 
normalized by a growth factor S and in our analysis we 
taken this into account. These data are not as useful as 
the redshift distortion factors for the determination of 7 
for two reasons 

1. The rms fluctuation ag(z) is not connected directly 
with the growth rate f(z). Instead, it is related 
with f(z) through the integral of eq. (|2. 12[) . 

2. Most of the Ly — a as data appear at high redshifts 
where ACDM is degenerate with most other dark 
energy models. 

Using the data of Table II we construct the correspond- 
ing xl defined as 



^ 0m ,7) = E 



Sobs(zi, Zi+l) — Sth(zi, 2i+l) 



(2.23) 

where <r Sobts 4 is derived by error propagation from the 
corresponding la errors of as(zi) and as(zi + \) while 
Sth(zi, Zi + i) is defined in eq. (|2.22|) . We can thus con- 
struct the combined Xtot(^0m,7) as 

xL(^o,„,7) = X/(fiom,7) +X«(fiom,7) (2-24) 



Notice however that all Refs report their data points 
having assumed a ACDM model when converting red- 
shifts to distances with fi 0m = 0.3 (except [27[ that used 
Q 0m = 0.25 and [H,[ll that used O 0m = 0.26). 

Setting f2om = 0.3 and minimizing Xtot w hh respect to 
7 we find 



7 = 0.674 



+0.195 
0.169 



(2.25) 



which differs somewhat from the corresponding result of 
Ref. because we have used a broader dataset, a dif- 
ferent parametrization for / and we have assumed ACDM 
as our fiducial model thus avoiding the marginalization 
of the parameter wq. The result (|2.25p indicates that the 
ACDM value of 7 = ^- = 0.545 is well within la from 
the best fit and is clearly consistent with data. The im- 
posed constraints however are rather weak and even a 
flat model with matter only (SCDM) predicting 7 = 0.6 
(set w = in eq. (|2.10j) ) is consistent with the data. 
Also, this result indicates that if it was only the value of 
fobs and the measured f3 that would be required then it 
could have been obtained trivially from the ACDM f Q b s 
(which is analytically known) and the measured (3 with- 
out need for a simulation. In that case we would also 
have found almost perfect agreement with ACDM. In- 
stead we find a somewhat larger value of 7. Performing 
the same analysis but by excluding the three datapoints 
at z = 0.35,0.55 and 1.4 yields a slightly different value 
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FIG. 2: The cosmological data for the growth rate f(z) 
along with the best theoretical fit / = f2 m (z) 7 with S7o m = 
0.3 (black continuous line) and the corresponding la errors 
(shaded region). The errorboxes on f are obtained using the 
ratios at the specific redshifts. Clearly, the best fit shows a 
minor difference from ACDM (blue dashed line) only at low 
redshifts. 



for 



7 = 0.73 ±0.23 



(2.26) 



but again well within la from the ACDM value. 

Ignoring the as(z) data of Table II leads to a negligible 
change in the best fit of all 22 points to 7 = 0.663 ± 0.2. 
This is consistent with the above discussion on the use- 
fulness of these data. Alternatively, assuming the ACDM 
value for 7 (7 = yj) and minimizing with respect to f2om 
we find Q 0m = 0.24±g;$. 

The cosmological data for the growth rate f(z) are 
shown in Fig. 2 along with the best theoretical fit / = 
Qm(z) 7 with f^om = 0.3 and the corresponding la errors 
(shaded region). In the same plot we show (dashed line) 
the corresponding Jacdm{z) obtained by solving numer- 
ically eq. (|2.4[) for w = — 1 and f^om = 0.3. Clearly, the 
best fit (continuous line) shows a minor difference from 
ACDM (dashed line) only at low redshifts. We therefore 
conclude that ACDM is consistent with current data for 
the growth factor 8{z) and the consistency is maximized 
for f2 0m 



24+ 09 



III. ALTERNATIVE TESTS OF ACDM USING 

«(*) 

The test of ACDM discussed in the previous section 
requires prior knowledge of the parameter f2om since the 
growth function data are not accurate enough to produce 
a simultaneous fit for both fio m and 7 with reasonably 
small errors. Thus it is useful to derive a test that de- 
pends only on the observables H(z) and 5(z). Such a 
consistency test has been recently discussed in Ref. [ToT ] 
and involves both HQs), 5(z) and their derivatives (see 



FIG. 3: The la range (shaded region) of the lhs of eq. (|3.6[) 
(black continuous line). Notice that by construction it is 
independent of the redshift z since we have assumed that 
the geometric part of eq. (|3.6|l (H(z)) behaves like ACDM 
(fiom ~ 0.3). The value corresponding to ACDM for both 
H(z) and 8(z) (dashed line) is clearly well within la from the 
best fit (continuous line). This is to be expected since the 
range of 7 in eq. (12.25[) includes the value 7 = tt ■ 



eq. (jl.6p ). Here we derive an improved version of this 
test that is independent of the derivative of H(z) and 
therefore it is less subject to observational errors. We 
start from eq. (|2.ip and change variables from t to a to 
obtain 



dH(a) 
da 
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2 3tto m H Q 5 
a 5 S' 



where ' = 4-. The solution to (|3~Tj) is [2£ 
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which may be expressed in redshift space as 
H 2 _ 3n 0m (l + z) 2 f°° S(z)S'(z) 



Hi 



S'(z) 2 

and setting z = we find 



l + Z 



dz 



S(z)S'(z) 
l + z 



dz 



(3.1) 



(3.2) 



(3.3) 



(3.4) 



In order to promote eq. (|3.4[) to a null test for ACDM 
we must express 57o m in terms of geometrical observables 
like HQs). 

In the context of ACDM we have 



n 



ACDM 1 



H(z) 
Hn 



1 



(l + z) 3 -l 



(3.5) 



leading to 9^ DM {z) = fi 0m when H(z) has the ACDM 



6 



form. By dividing f|3. 5|) with (|3 .4|) we have 

o ( H(z) 2 ,\ roc S(z)S'(z) , 

nj CDM ( Z ) 1= H^r^H 1=n 

O 0m [(l±z) 3 -l]<5'(0) 2 

(3.6) 

where the last equality should hold only if ACDM is a 
valid theory. Using ([232]) . (|2~T3"j) and $tl5§ in we 
can find the la range of the lhs of eq. (|3.6p which is 
shown in Fig. 3 (shaded region). Notice that by con- 
struction it is independent of the redshift z since we have 
assumed that the geometric part of eq. (|3.6[) (H(z)) be- 
haves like ACDM (fiom = 0.3). The value correspond- 
ing to ACDM is clearly well within la as shown in Fig. 3. 
This is to be expected since the range of 7 in eq. (|2.25|) 
includes the value 7 = jr. Had we allowed for a more 
general form of H (z) then the la range of eq. (|3.6p would 
be redshift dependent. However, as discussed in section 
II, in that case the allowed range of 7 would be less re- 
liable since the surveys leading to the data of Tables I 
and II convert redshifts to distances using ACDM with 
fiom = 0.3. Therefore, these datapoints can only be used 
to test ACDM . 



easily tested by using geometrical tests alone. However, 
in that case, care should be taken to reanalyze the power 
spectra using the proper form of H (z) when converting 
redshifts to distances. 

Nevertheless, given the current uncertainties in the 
growth rate observations, data of much better quality 
will be needed in order to distinguish between ACDM 
and modified gravity theories. These data will most likely 
come from large scale weak lensing surveys like DUNE, 
which is expected to measure the the equation of state of 
dark energy to a precision better than 5% [301 ] . 

Numerical Analysis: The mathematica files with 
the numerical analysis of this stu dy may be found 
at |http://nesseris.physics.uoi.gr/growth/growth.htm| or 
may be sent by e-mail upon request. 
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V. APPENDIX 



IV. CONCLUSION 

We have compiled a dataset consisting of redshift dis- 
tortion factors (3{z) and rms mass fluctuations ag(z) at 
various redshifts obtained from galaxy and Ly — a for- 
est redshift surveys. Using this dataset we have obtained 
the best fit form of the linear growth function 5{z) us- 
ing the Wang and Steinhardt parametrization (|2.13p with 
flom = 0.3. We have found a best fit value of 7 as 
7 = 0.674to;!69, a range which includes the ACDM value 
7 = Jr. Thus ACDM is in excellent agreement with 
current linear growth data. We have reached the same 
conclusion by applying a generalized version of the null 
test of Ref . [l(| • 

The combination of geometrical (H(z)) and dynamical 
(<5(z)) observables used here to test ACDM could also be 
used to test modified gravity theories which can not be 



It is straightforward to generalize our analysis of sec- 
tion II using (|2~T4"j) in (|2~T3"|) . ([237]) . ([2~2"2"|) . In that case 
we have two parameters to fit and we obtain 

70 = 0.77 ±0.29 (5.1) 
7 = -0.38 ±0.85 (5.2) 

Clearly, the la error region is significantly increased in 
that case (see Fig. 4) due to the introduction of the addi- 
tional parameter and ACDM (70 = 7 ~ 0) remains 
well within la from the best fit. 

In order to avoid the increased error region we may 
also set 70 = ji and fit 7 only. In this case we get 

7 = 0.17 ±0.54 (5.3) 

which is again consistent with ACDM with more reason- 
able errors. 
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